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1. INTRODUCTION 


Surface catalyzed reactions occurring in open systems exhibit various ki- 
netic (or nonequilibrium) phase transitions' *. In general such processes are very 
complex and strongly affected by adsorption, desorption, diffusion and reaction 
rates^^. One of the most intensively studied systems is the lattice model for the 
carbon monoxide oxidation on a single-crystal catalyst surface" ^. Experi- 
mental studies'*'* have demonstrated that the carbon monoxide oxidation is a 
Langmuir-Hinshelwood process, and its basic steps are: i) adsorption of CO 
molecule on a single vacant site, 11) dissociative adsorption of O; molecule on a 
pair of adjacent empty sites, and iii) reaction between CO and O followed by 
desorption of the CO; molecule. Since the binding energy for CO is lower than 
the binding energy of oxygen atoms'””°, the former undergoes surface diffusion 
during the process, while oxygen atoms are considerably less mobile. 

From the numerous mean-field**'*'** and Monte Carlo**77 studies it 
follows that this system exhibits two kinetic phase transitions as the impinge- 
ment rates of the reactants change. The model introduced by Ziff et al. assumes 
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the infinite reaction rate and no diffusion of the reactants. These authors have 
performed extensive Monte Carlo simulation for the impingement rate of CO 
molecules pco between O and 1 (the corresponding impingement rate for O» 
molecules is equal to pg, =1— Poo - At low pco the poisoning of the surface by 


oxygen atoms occurs. As the probability of CO adsorption increases, however, 
the continuous transition between this O — poisoned and the reactive states oc- 
curs at pt, equal to 0.389. The reaction window spans the region of pco be- 


tween poo and pe, =0.525. At pco equal to po the second, first-order transi- 


tion occurs and for higher values of pco the CO-poisoned state appears. The 
locations of these phase transitions change when the diffusion of reactants is 
allowed? ^6, when the reaction rates are finite? as well as when the desorption 
of reactants can take place*”'. 

Another interesting problem is the influence of the catalyst surface nonuni- 
formity on the reactions and the kinetic phase transitions. The presence of sur- 
face defects on the properties of reacting systems has been previously consid- 
ered by many authors”. In particular, Monte Carlo simulations performed on 
percolating clusters and Sierpinsky carpet”***', both with the fractal dimension 
Dpr = 1.89 have shown that the monomer-dimer reaction system exhibits two 
second-order (continuous) transitions, instead of one first-order and one second- 
-order transition found for the two-dimensional systems (with Dp = 2). Effects 
of surface deactivation on the catalytic reactions have also been studied for the 
oxidation of methanol to CO on sulfur containing Fe(100)??* and Ni(100) * 
surfaces. It has been shown that as the surface coverage of sulfur increases the 
production of CO considerably decreases and another reaction, leading to the 
formation of formaldehyde, takes precedence. 

In this paper we consider a simple monomer-dimer reaction and discuss the 
effects of the reaction rate magnitude on the locations of the kinetic phase tran- 
sitions and the changes in the mechanism of these transitions when the surface 
contains randomly distributed nonactive impurities. 


2. THE MODEL AND MONTE CARLO METHOD 


We consider a simple square lattice of the size L x L and assume that each 
lattice site is accessible for adsorbing molecules. The adsorbing mixture consists 
of two types of molecules A> (O5) and B (CO). In general, we assume that only a 
certain fraction y, of the total number of sites (2) is active. These active sites 
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are assumed to be distributed randomly over the lattice. The particles B arrive at 
the surface with the probability pg and the particles A; with the probability 
Pa 71— p, . Adsorption of B occurs whenever the particle meets an empty site, 


while adsorption of A; requires a pair of adjacent empty sites. The state of the 
adsorbed B molecules is not affected by the type of adsorption site (active — 
marked by S$”, and nonactive — marked by S). On the other hand, the state of the 
adsorbed A; molecule depends on the type of sites. If both adsorbing sites are 
nonactive the A; molecule does not dissociate. When at least one of the adsorb- 
ing sites is active, however, the adsorption of A; is accompanied by dissociation. 
Only these atoms A which are adsorbed on active sites are able to react with B 
particles and form the product AB (CO;) which immediately desorbs from the 
surface leaving two empty sites. The reaction rate is assumed to be finite and 
equal to k. Particles B, are assumed to be able to diffuse over the surface at the 
jump rate [p independently of the lattice site type. Also the molecules A> and 
the atoms A adsorbed on nonactive sites can travel over the surface at the jump 
rates equal to T, and T4, respectively. The atoms A adsorbed on active sites 


(marked by A;) remain localized, however. Whenever the diffusing adsorbed 


molecule A; meets an active site it dissociates into two A atoms. 
Thus, the adsorption processes can be represented as 


B! + S—> B" (1) 
A? +25—> B; (2) 

when both sites are nonactive, 
As +25 ——>2A” (3) 


when both adsorbing sites are active (A" denotes adsorbed activated A atoms) 
and 


A: $4 S' —A' e A" (4) 


when only one of the adsorbing sites is active. The dissociative diffusion of A» 
molecules can be represented as 
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A; +S — A c A" +S (5) 
Finally, the reaction occurs according to the equation 
A" 4- B* —3À5 AB* 4 S" +S(S”) (6) 


and the product immediately leaves the surface. 

To study the behavior of such systems we have used the Monte Carlo 
method. The applied algorithm is quite similar to that used previously by Sholl 
and Skodje. Thus, we assume that a single Monte Carlo attempt can be adsorp- 
tion of either A» or B particle (with the probabilities given by p, and 


Ps 71- p, , respectively, diffusion of B, diffusion of nondissocited A; mole- 


cule or A atom adsorbed on nonactive sites, or the reaction between A and B, 
accompanied by the desorption of the resulting AB molecule. Thus, a single 
Monte Carlo step defines a time unit as equal to 


t, =1+T, +I +I, +2k (7) 


Monte Carlo attempts. Calculations have been carried out for lattices of the size 
L= 120 and for y, ranging between 0.025 and 1. The number of Monte Carlo 
steps used in this study varied considerably with the parameters describing the 
process. In general, B poisoning occurs much more rapidly than the processes 
occurring in the reaction window and A poisoning’. In the case of reactive states 
the simulation required up to 10° moves before the steady state was reached. 
Then, another 10* Monte Carlo steps has been performed to calculate averages. 
In all cases, several independent runs have been performed and the final results 
have been obtained by averaging over these runs. In the case of partially nonac- 
tive surfaces, the surface configuration of active sites has been generated sepa- 
rately in each run. Standard periodic boundary conditions have been applied in 
both directions. 


3. THE EFFECT OF REACTION RATE ON THE PHASE DIAGRAM 


Here we consider the effect of the changes in the reaction rate on the loca- 
tion of both kinetic phase transitions in the case of completely active surface 
(Ya = 1). Although the problems of difussion and reaction rates have been al- 
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ready considered in the literature" ****"*, we have decided to include a brief 


discussion for completeness. The performed calculations have also served as a 
check for the correctness of our simulation algorithm and code. 

In the systems with fully active surface, all the adsorbed A; molecules dis- 
sociate and are completely localized, while the molecules B are assumed to 
travel over the surface with the jump rate [p = 1. Fig. 1 shows the changes in the 
reaction window, i.e., the locations of the kinetic phase transitions (Pó, and 


Pi, , versus the reaction rate k. The result for k = œ has been taken from Ref. (4) 


and corresponds to the situation in which both reactants are fully localized. 


—— u " 
0.2 0.4 0.6 0.8 
1/k 


Fig. 1. Dependence of the kinetic phase transition points on the reciprocal reaction rate (1/k) 
for the systems with [p = 1. Only the results for k = o correspond to I’; = 0 
and have been taken from Ref. (4) 
Zależność między punktami kinetycznej przemiany fazowej i odwrotnością 
stałej szybkości reakcji (1/k) w przypadku układów z I5 = 1. Jedynie wyniki 
otrzymane dla k = © odpowiadają I's = 0 i zostały zaczerpnięte z pracy (4) 


We observe that as the reaction rate becomes smaller the reaction window 
becomes narrower, mostly due to rather large change in the location of the first 
order kinetic phase transition Pu . The continuous transition separating reactive 


and B-poisoned states is only very slightly affected by the reaction rate. Also, 
the diffusion of reactants has a small influence on the reaction window when the 
reaction rate is high. On the other hand, for small reaction rates the increase in 
the jump rate of the reacting particles is expected to be more important". 

In Fig. 2 we present the examples of phase diagrams for the systems char- 
acterized by k=1, 2 and 4, and the jump rate for the component B equal to 


Fs =1. Our results for k=1 are the same as reported earlier by Sholl and 
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Fig. 2. Phase diagrams for the systems with k = 1 (part (a)), k = 2 (part (b)) 
and k = 4 (part (c)), obtained for the model with I’; = 1 
Diagramy fazowe odpowiadające układom z k= 1 (część a), k = 2 (część b) 
oraz k = 4 (część c) odpowiadające modelowi z Tg = 1 
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Skodje®. In particular, we have found that the kinetic data in the region of B- 
-poisoning transition approximately obey the scaling relations 


e,()- Flo? - pa x] (8) 
e,()- Flo? - p, v] (9) 


where PA, represents the location of the spinodal point. 


The results for the production rate, presented in Fig. 2, require some com- 
ments. The production rates have been expressed per single Monte Carlo move 
and not per single Monte Carlo step, consisting of t, MC moves. As the reaction 
rate constant k increases it leads to the increase of t,. Thus, the number of reac- 
tion events per time unit changes with k. In order to express the production rates 
per MC step, the results must be multiplied by ¢,. Fig. 3 presents the production 
rates just at the onset of the reaction window (at P^, ) for systems with different 


k and expressed per Monte Carlo step. Now, we observe that the production rate 
increases with k as expected, going smoothly to the value obtained for k = oo. 


0.25 F> 
PRA 


0.20 


0.15 


Fig. 3. The maximum production rate at the onset of the reaction window at Pr 


as a function of the reciprocal reaction rate, for the systems with Tg = 1. 
Results for k = o have been taken from Ref. (4) 
Maksima szybkości produkcji w punktach początkowych okna reakcji P^ 


jako funkcja odwrotności stałej szybkości reakcji w przypadku układów z Tg= 1. 
Wyniki odpowiadające k = w zaczerpnięto z pracy (4) 
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4. REACTION ON PARTIALLY DEACTIVATED SURFACE 


For small values of p, (lower than pj ) and when all sites are active the 


surface becomes B-poisoned, i.e., entirely covered by the adsorbed B particles. 
When the fraction of active sites becomes lower than unity, we also observe 
poisoning, but with nonzero concentrations of B and both A> and A particles. 
This is illustrated in Fig. 4, which presents the surface coverages of A, A; and B 
particles in the steady state for surfaces characterized by different fraction of 
nonactive impurity sites (Yna = 1 — Yna). The total surface coverage is very close 
to unity in all cases, however. Note that the surface coverage of nondissociated 
A» molecules is given by the ratio 6, =N, / D , where N 4, 1$ the number of A> 


molecules in the system. Thus, the number of sites occupied by A» molecules is 
equal to 2N, . The presence of A atoms results from the fact that some of them 


are trapped on nonactive sites and surrounded by either nondissociated A; or B 
molecules. According to the assumed model, the atoms A adsorbed on nonactive 
sites do not react with B molecules. As the fraction of nonactive sites goes to 
unity, the reaction ceases completely and the surface becomes covered by A; and 
B particles in the proportions fully determined by p, . It is obvious that the 


considered here systems with partially deactivated surface are examples of the 
well-known models exhibiting multiple absorbing states ™*®, 

In the case of fully active surface the model predicts two kinetic phase tran- 
sitions as p, changes. One is the first-order transition separating the B- 


-poisoned and reactive states and it is located at fs; and the second is the con- 


tinuous transition separating the reactive and A-poisoned states and located at 
DA, . When the nonactive impurities are present in the system, we find that both 


transitions become continuous. Fig. 5 presents examples of phase diagrams for 
the systems with different fractions of nonactive impurities. An interesting fea- 
ture of partially deactivated systems is a considerable widening of the reaction 
window, though the production rate of AB molecules becomes lower and lower 
as the fraction of nonactive sites increases. When the production rate is related 
to the number of active sites only and not to the total number of sites, we find 
that the maxima of the production rate for partially deactivated surfaces appear 
at different values of p, for a given fraction of active sites y, 
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Fig. 4. Surface coverages of species A (part (a)), B (part (b)) and A> (part (c)) 
for partially deactivated surfaces and different values of p, below the first-order 


kinetic phase transition, for the systems with k = 1, T4 = Tg =T, = 1 


Stopnie pokrycia powierzchni substancji A (część a), B (część b) oraz A> (część c) odpowiadające 
częściowo deaktywowanym powierzchniom i różnym wartościom p, poniżej punktu kinetycznej 


przemiany fazowej pierwszego rzędu w przypadku układów z k = 1, 4-2 sz T4, = 1 


ANDRZEJ PATRYKIEJEW, STEFAN SOKOŁOWSKI 


Pr 


0.02 


0.01 


Pr 


0.02 


0.01 


Pr 


0.02 


0.01 


Fig. 5. Examples of phase diagrams for systems with partially deactivated surface 
with y, = 0.9 (a), 0.7 (b) and 0.5 (c). All results have been obtained 
for k= 1 andTq=Tp=l, =1 
Przykłady diagramów fazowych odpowiadających układom z częściowo 
deaktywowaną powierzchnią z y, = 0,9 (a), 0,7 (b) i 0,5 (c). Wszystkie wyniki uzyskano 
w przypadku k= Li 42 Fg zI, =1 
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Fig. 6 presents the example of the results obtained for p, equal to 0.6. 
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Fig. 6. The changes in the surface coverages of different species versus the degree of surface deac- 
tivation for p, = 0.6. The inset shows the production rate (per Monte Carlo move) 
Zmiany stopnia pokrycia powierzchni róznych substancji wzgledem stopnia deaktywacji 
powierzchni w przypadku p, = 0.6. Wstawka ukazuje szybkości produkcji 
(na jeden krok Monte Carlo) 


We observe that initially, for Yn between zero and about 0.225, surface deacti- 
vation leads to nearly linear increase of the production rate and also causes 
nearly linear changes of the surface coverages for different species. Then, as Yna 
increases, the production rate gradually decreases to zero, and the surface be- 
comes poisoned for Yna ~ 0.42. Subsequent changes of the coverages for all spe- 
cies are quite similar to the results given in Fig. 4. From the plot of the maxi- 
mum production rate versus p, (see Fig. 7) we find that the most optimum 


regime for the production of the reaction product AB corresponds to p,, = 0.65 


and y, = 0.8. 
In general, surface deactivation shifts the locations of both transitions, and 
hence the reaction window, towards higher values of p, (see Fig. 8). We find 


that the reaction window becomes wider for partially deactivated surface, 
though the AB production rate rapidly decreases when the degree of surface 
deactivation increases. 
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Fig. 7. The maximum production rate for partially deactivated systems 
as a function of p, (k=l,T4=[p=1, 7 1) 
Maksima szybkości produkcji w przypadku częściowo deaktywowanych powierzchni 
jako funkcja p, (k=l,T4=[5=T, = 1) 
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Fig. 8. Locations of the transition points in partially deactivated systems 
with k= 1 andl4=T5=T, =1. 


0 and o represent the results obtained for k= 4andTq=Tp=T, =1 


Położenie punktów przemiany fazowej w przypadku układów z częściowo 
deaktywowaną powierzchnią przy k= 1iT4=Tp=T, =1. 


0 i o odpowiadają wynikom uzyskanym w przypadku k=4iI4=Ip=T, =1 
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An interesting problem that has not been considered here in detail, is the 
crossover between the first-order and continuous transition as the concentration 
of nonactive impurities in the system increases. From theoretical studies of 
equilibrium phase transitions in lattice systems with impurities^^' it follows 
that two scenarios are possible. For example, in the case of the random field 
Ising model one expects that even infinitesimal randomness destroys any sharp 
phase transition??. On the other hand, computer simulations for the adsorbed 
films (also modeled by the lattice gas model) formed on lattices containing two 
types of sites have shown that for sufficiently small concentration of less (or 
more) active impurity sites, the transition is still present^'. The situation consid- 
ered in this work is more similar to this model than to the above mentioned ran- 
dom field Ising model. Therefore, one can expect that the crossover between the 
first-order and continuous transition is likely to occur only for sufficiently high 
concentration of nonactive sites. Our results obtained for Yn = 0.05 do not show 
any rounding of the first-order kinetic transition, while for Yna — 0.1 the transi- 
tion appears to be continuous. 

Another emerging question is the effect of impurities on the critical behav- 
ior of the continuous kinetic phase transition occurring at Ph . In the case of 


fully active surface and infinite reaction rate, the critical properties of the con- 
sidered model are well known?™™. It has been found that this transition belongs 
to the universality class of the directed percolation model”. The order parameter 
for this transition, 05, is expected to behave as 


8,(p,)-« lp, - pil qo) 


with the exponent B equal to 0.58. One can expect that surface deactivation 
should considerably affect this transition, or destroy it completely. Our results 
suggest that for sufficiently highly deactivated surface this transition becomes 
gradual, i.e., it is not any true phase transition any more (cf Fig. 5). In particular, 
it seems that the power-law given by eqn. (10) is not obeyed. In order to address 
this problem, and obtain more reliable results, the simulations for much larger 
systems are needed, however. In particular, the method developed to study dam- 
age spreading at equilibrium phase transitions", seems quite well suited for 
such problems? and will be considered in our future work, where we shall also 
present the solutions of the mean field equations for the model considered here. 
Although, the considered model is extremely simple and the presented re- 
sults have been obtained for rather small systems, but we believe that it is a 


68 ANDRZEJ PATRYKIEJEW, STEFAN SOKOŁOWSKI 


good starting point for investigation of more complex systems. In particular, it is 
possible to consider the situations in which there are several possible reactions 
occurring at the surface which compete one with another depending on the state 
of adsorbed species. Our preliminary results demonstrate that the production 
rate strongly depends on the surface activity and that surface deactivation may 
lead to an increase in the reaction yield (per one active site) in some cases. It is 
probably due to the assumed mobility of dimer molecules adsorbed on nonactive 
sites. The considered situation is somewhat similar to the model allowing diffu- 
sion of the dissociated A; particles over the surface. We shall return to this 
problem in our future work. 
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STRESZCZENIE 


Przedstawiono wyniki symulacji Monte Carlo dla modelowej, katalizowanej powierzchniowo 
reakcji A; + 2B — 2AB. Rozważono wpływ zmian szykości reakcji oraz deaktywacji powierzchni 
na kinetyczne przemiany fazowe. Pokazano, że deaktywacja powierzchni katalizatora w istotny 
sposób wpływa na położenie okna reakcji oraz na mechanizm obserwowanych przemian 
fazowych. W szczególności, pokazano, że deaktywacja powierzchni zmienia rząd kinetycznej 
przemiany fazowej w kierunku przemiany ciagłej. 


